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ABSTRACT 

The Carina Nebula is one of the youngest, most active sites of massive star formation in our 
Galaxy. In this nebula, we have discovered a bright X-ray source that has persisted for ~30 years. 
The soft X-ray spectrum, consistent with kT ~128 eV blackbody radiation with mild extinction, 
and no counterpart in the near- and mid-infrared wavelengths indicate that it is a ~10^-year-old 
neutron star housed in the Carina Nebula. Current star formation theory does not suggest that 
the progenitor of the neutron star and massive stars in the Carina Nebula, in particular rj Car, 
are coeval. This result suggests that the Carina Nebula experienced at least two major episodes 
of massive star formation. The neutron star may be responsible for remnants of high energy 
activity seen in multiple wavelengths. 



Subject headings: supernova remnants 
bubbles — X-rays: stars 

1. Introduction 



Massive stars (M > IOMq) are born from gi- 
ant molecular clouds along with many lower mass 
stars, forming a stellar cluster or association. Mas- 
sive stars evolve orders of magnitude more quickly 
than lower mass stars, and die through super- 
nova or hypernova explosions in ^10^ years. High 
interstellar pressures that massive stars produce 
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through their strong UV radiation and supernova 
explosions can compress a pre-existing cloud and 
trigg er the formation of new stars (cf. lElmegreen 



19981) 



The Carina Nebula is one of the most massive 
star forming regions in our Galaxy. It contains 
two massive stellar clusters, Trumpler 14 and 16 
(Figure 1), possessing over 50 massive sta rs with 



spect ral types earher than 06 (>40 M©, I Smith 
l2006h . The nebula is also home to one of the most 
massive stars in our Galaxy, rj Carinae, which has 
an estimated initial mas s ^,150 M© and a current 



mass of about 90 M© ( Hilher et al.ll200lh 



At a distance of only ~2.3 kpc ( Davidson fc Humphreys! 
19971) . the Carina Nebula is one of the best sites 



for studying how very massive stars form and af- 
fect their environment. It shows signatures of 
violent activities: a bipolar supershell structure 



([Smith et al.ll2000r). strong turbul ence in interstel- 
lar clouds ( Yonekura et al. 2005f ). a nd hot X-ray 



plasma through t he en ti re nebula (ISeward et al 



19791: iTownslevI 120061 : iHamaguchi et al 



20071 : 



Ezoe et ahf 2008 ) . Two primary mechanisms have 



been proposed to produce these structures. One is 
strong winds and UV radiation from massive stars 
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in the nebula, while the other is supernova ex- 
plosions. The energy budget and elemental abun- 
dance distribution favo r the supernova mechanism 



( Yonekura et al. l2005t iHamaguchi et all l2007t ) 



However, neither black hole, neutron star, nor 
clear remnant from a supernova, has been found 
in the Carina Nebula. In addition, any massive 
supernova progenitor is unlikely to be co-eval with 
the observed stars since it would be more massive 
than 7] Car; this means that, if a supernova oc- 
curred in the Carina Nebula, its progenitor would 
probably have formed much earlier than 77 Car. 
Thus, detection of any compact object in the Ca- 
rina Nebula would be of vital interest for under- 
standing the star formation history of the region. 

From multiple sets of X-ray data, we found 
a promising neutron star candidate at the heart 
of the Cari na Nebula. The source, which we 
call EHG7 (|Ezoe et all I2OO8I ). was detected as 



a very soft source in X-ray images 8.5' south- 
east of 77 Car, equivalent to 5.7 pc in the pro- 
jected physical distance assuming d =2.3 kpc (See 
Figure 1). The X-ray source has also been re- 
ported contemporaneous l y in f our other papers 
(lAlbacete Colombo et a"l]|2003l: lEzoe et al.| [2008: 



Pires fc MotchI l2008t IPires et all I2OO8I ). the last 
two also suggested EHG7 as a neutron star based 
on XMM-Newton and Chandra data. This paper 
presents a comprehensive study of the character- 
istics of EHG7 using all available X-ray data from 
Einstein and ROSAT observations to the latest 
XMM-Newton observations in 2009 to present con- 
clusive evidence of the compact nature of EHG7. 

2. X-ray Observations and Analysis 

We found 32 X-ray observations with EHG7 in 
the field of view and spatial resolution better than 
1' in the HEASARC archive and our proprietary 
data (Table [T|) . Observations with lower spatial 
resolution (e.g. Einstein IPC, ASCA) failed to de- 
tect EHG7 due to severe contamination from sur- 
rounding soft diffuse X-ray emission. Throughout 
this document, individual observations with Ein- 
stein, ROSAT, XMM-Newton, and Chandra obsei- 
vations are designated EIN, ROS, XMM, CXO re- 
spectively, subscripted with the year, month, and 
day of the start time of the observations. 

We used the software package SAS, version 
8.0.1 for analysis of the XMM-Newton EPIC data 



(jStriider et al.l 120011 ). We followed the standard 
methocfl hi data processing — screening out high 
background periods, selecting photon events and 
generating spectral response files. For the tim- 
ing and spectral analysis, we took source events 
from an encircled region with 30" radius and back- 
ground events from an annulus region with 60" 
outer and 30" inner radii. For data sets with the 
source close to 2 CCD gaps (see Table [J), source 
events are collected only from a CCD chip with 
the source peak. Comparison of these spectra with 
spectra of the other cameras in the same observa- 
tion looked consistent, but we did not use them 
for flux measurements. 

For source position determination on the Chan- 
dra image, we used a custom anal ysis method 
for C handra data sets developed by iBroos et al. 
(|2007l ) and Broos et al. (in prep). In producing 



the Chandra spectral and timing data sets, we ana- 
lyzed level 2 event data using the software package 
CI AO, version 4.0. We took source events from an 
encircled region with a 8" radius and background 
events from an annulus region with 16" outer and 
8" inner radii. 

For general light curve and spectral analysis 
of the above data sets and for the ROSAT and 
Einstein data processing, we used the software 
package HEAsoft version 6.5.1. We started from 
the pipeline products of the ROSAT and Einstein 
data. For the ROSATima.gc analysis, we excluded 
PHA channels above 10 to reduce particle back- 
ground events for the HRI and extracted events 
between 0.3-2 keV for the PSPC. For the PSPC 
spectral and timing analysis, we took source events 
from an encircled region with 30" radius and back- 
ground events from an annulus region with 60" 
outer and 30" inner radii. For Einstein HRI cal- 
ibration information, wc referred the mirror vi- 
gnetting at the location of EHG7 in EIN781224 (8.5' 
off-axis) to lHarnden et all (|l984l ) (91% at 0.6 keV) 
and estimated encircled energy of the point spread 
function (PSF) within a 32"x32" box at 76% us- 
ing a calibration observation of Cyg X-1 at a simi- 
lar off-axis angle (10': the dataset H1956N35.xic). 
The dead time was negligible and therefore ig- 
nored. 



^http;/ /xmm.esac.esa.int/sas/8.0.0/documentation/threads/ 
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Table 1 

X-RAY Observations with EHG7 in the Field of View 



Abbreviation 


Observation ID 


Observation Start 


Detector/Filter 


Offset 
(arcmin) 


Exposure 
(kscc) 


Flux (0.3-2 keV) 
(10^^^ orgs cm"^ s"^) 


EIN781214 


1074 


1978 


Dec. 14, 05:44 


HRI-2 


8.5 


17.4 


0.89±0.54 


ROSg00727 


CA150037H.N1 


1990 


Jul. z7, 0U:12 


HRI 


15.0 


3.3 


<1.1 


ROS911215 


US200108P.N1 


1991 


Dec. 15, 9:58 


PSPCB 


19.4 


1.6 


<2.2 


ROS92O6I2 


US900176P.N1 


1992 


Jun. 12, 22:33 


PSPCB 


8.7 


23.6 


0.94±0.16 


KUO92073I 


VVVjyUUooOrl.iNl 


1992 


Jul. 31, 01:03 


rlrti 


Q 7 
O. ( 


11/1 


1 r\A -1-/1 T ^ 


ROS92O8O9 


WG201262P.N1 


1992 


Aug. 9, 22:21 


PSPCB 


13.8 


5.7 


1.43±0.47 


ROS92O8IO 


WG200709P-1.N1 


1992 


Aug. 10, 19:03 


PSPCB 


35.4 


5.9 


<2.4 


ROS921215 


US900176P-1.N1 


1992 


Dec. 15, 17:39 


PSPCB 


8.7 


14.1 


0.69±0.18 


ROS940106 


WG900385H-2.N1 


1994 


Jan 6, 1:54 


HRI 


8.7 


0.5 


<7.7 


ROS94072I 


WG900385H-3.N1 


1994 


Jul. 21, 02:04 


HRI 


8.7 


40.1 


0.91±0.18 


ROS96O8I3 


US900644H.N1 


1996 


Aug. 13, 21:07 


HRI 


8.7 


1.7 


<3.3 


ROS971223 


US202331H.N1 


1997 


Dec. 23, 08:59 


HRI 


8.7 


46.5 


0.82±0.16 


XMMoo0726 


112580601 


2000 


Jul. 26, 05:08 


MOS/thick 


8.7 


-/33. 2/30. it 


1.04±0.18 


XMMoo0727 


112580701 


2000 


Jul. 27, 23:58 


MOS/thick 


8.7 


-/10.9/7.9t 


1.21±0.30 


XMM030125 


145740101 


2003 


Jan. 25, 12:58 


MOS/thick 


8.4 


-/6.9/6.9 + 


1.40±0.36 


XMMo30127A 


145740201 


2003 


Jan. 27, 01:04 


MOS/thick 


8.4 


-/6.8V6.9+ 


(1.21±0.25) 


XMMo30127B 


145740301 


2003 


Jan. 27, 20:37 


MOS/thick 


8.4 


-/6.8V6.8+ 


(1.11±0.23) 


XMMo30129A 


145740401 


2003 


Jan. 29, 01:41 


MOS/thick 


8.4 


-/8.4t/8.4 


0.95±0.27 


XMM03OI29B 


145740501 


2003 


Jan. 29, 23:55 


MOS/thick 


8.4 


-/6.9V6.9 


1.30±0.32 


XMM03O6O8 


160160101 


2003 


Jun. 08, 13:30 


MOS/thick 


8.4 


-/17.0/16.5 


1.10±0.19 


XMM03O6I3 


160160901 


2003 


Jun. 13, 23:52 


MOS/thick 


8.4 


-/31. 1/31.1 


1.03±0.15 


XMMo30722 


145780101 


2003 


Jul. 22, 01:51 


MOS/thick 


8.4 


-/8.4/8.4+ 


1.19±0.34 


XMM03O8O2 


160560101 


2003 


Aug. 02, 21:01 


MOS2/thick 


8.4 


-/-/11.9 


1.14±0.26 


XMM03O8O9 


160560201 


2003 


Aug. 09, 01:44 


MOS/thick 


8.4 


-/12.2/12.2 


1.06±0.20 


XMM03O8I8 


160560301 


2003 


Aug. 18, 15:23 


MOS/thick 


8.4 


-/18.5/18.5 


1.06±0.17 


XMMo41207 


206010101 


2004 


Dec. 7, 7:30 


pn/mcdium 


16.3 


19.4V-/- 


(0.88±0.14) 


XMM06OI3I 


311990101 


2006 


Jan. 31, 18:04 


pn&MOS/thick 


8.7 


25.2/65.3V65.3 


1.05±0.12 


CXO08O9O5 


9488 


2008 


Sep. 5, 21:24 


ACIS-I 


6.7 


59.4 


0.87±0.08 


XMM09OIO5 


560580101 


2009 


Jan 5, 10:23 


MOS/thick 


8.5 


-/14.0/13.9 


1.17±0.19 


XMM09OIO9 


560580201 


2009 


Jan 9, 14:28 


MOS/thick 


8.5 


-/II. 4/11. 4 


1.17±0.21 


XMM09OII5 


560580301 


2009 


Jan 15, 11:23 


MOS/thick 


8.5 


-/25.4/25.4 


0.97±0.15 


XMM090202 


560580401 


2009 


Feb 2, 04:46 


MOS2/thick 


8.5 


-/-/26.2 


1.09±0.18 



Note.— Abbreviation: EIN: Einstein, ROS: ROSAT, XMM: XMM-Newton, CXO: Chandra, Detector/Filter: optical blocking 
filter for XMM-Newton, Offset: angle of EHG7 from the nominal coordinate, Exposure: pn/MOSl/MOS2 for XMM-Newton, 
Flux: Error denotes the 90% confidence range. The upper limit is at a 3(7 level. The fiux errors include uncertainty of 
the absolute fiux calibration of each telescope (Einstein, ROSAT, XMM-Newton: 10%, Chandra: 5%, ROSAT User's Hand- 
book — http://heasarc.gsfc.nasa.gov/docs/rosat/ruh/handbook/handbook.html [Beuermann et al.| ( |2006f) , XMM-SOC-CAL-TN-0018 — 
[http:/ /xmm2.esac.esa.int/external/xmm_sw_cal/calib/documentation/index.shtml9^EPIC Calibration Requirements and Present Status 
— |http:/ / asc. harvard. edu/cal/docs/cal_present_status.html:j^abs_eff| l. ^the source peak was close to 2 CCD gaps or an edge of the field of 
view, 'i'hose data sets were not used for fiux measurements in the right column. 



3 




Fig. 1. — The central part of the Carma Nebula in X-rays taken with XMM-Newton {Left, red: 0.4—0.7 keV, 
green: 0.7—1.3 keV, blue: 2—7 keV) and in optical emission lines {Right, Credit: Nathan Smith, University 
of Minnesota/NOAO/AURA/NSF). EHG7, the source discussed here, is denoted by the white arrow. It 
is located on the conspicuous V-shaped dark lane running across the optical image and is surrounded by 
diffuse X-ray emission. The black dotted line and circles depict positions of the super massive star rj Car 
and massive stellar clusters, Trumplcr 14 and 16, respectively. 



3. Results 

The Chandra observation CXOo80905: per- 
formed as a part of a Very Large Project to map 
the Carina Nebula in X-rays (PI: Leisa Towns- 
ley), gives the most precise and reliable abso- 
lute position of EIIG7 (Figure d]). In this ob- 
servation, EIIG7 was imaged 6.7' off-axis on the 
ACIS-I detector. The absolute coordinates of the 
source were determined to be (a20oo, ^2000) = 
(10'>46'°08.72", -59°43'06.5") with less than 1" 
positional uncertainty after correcting the Chan- 
dra astrometry by cross correlation of X-ray and 
infrared sources from the 2MASS catalogue and 
correlating the PSF at the location of this source 
with the image of the source. 

We searched for near- and mid-infrared coun- 
terparts from images obtained with the IRSF 
(InfraRed Survey Facili ty) 1.4m telescopes (see 



Sanchawala et al. 2007allbl . for details) and in all 



available Svitzer I nfraRed Array Camera (IRAC; 



Fazio et al.l 120041 ) observations combined, but 



the error circle of the X-ray source (Figure [2|). We 
used standard aperture photometry routines in 
IRAF to estimate the S-tr upper-limits, which were 
18.5 mag, 19.5 mag and 18.5 mag in the J,H and 
Ks bands and 0.84, 0.51, 4.7, and 13.6 mJy in the 
IRAC 3.6, 4.5, 5.8, and 8.0 /xm bands, respectively. 
Assuming blackbody spectra with temperatures 
between 2000 and 10000 K, the bolometric flux of 
EIIG7 was limited to <5xl0~^'^ ergs cm~^ s~^ 



these images showed no counterpart object within 



We note that Pires et al. ( 2008[ ) measured upper 
limits of 25™ in R and i?-band optical images. 

X-ray emission from EIIG7 was detected at 
above 3a significance in all observations in Table [T] 
except for 5 iJOS'^T' observations with short expo- 
sures. None of these observations showed any sig- 
nificant time variation from EIIG7 at greater than 
90% confidence in a fit of each light curve in the 
whole instrumental energy band with a constant 
flux model. For data sets with no spectral resolu- 
tion and/or poor photon statistics, we measured 
the photon count rate or upper limit in whole band 
images using the sosta package in ximage, taking 
a suitable nearby source free region as background. 
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10:46:10.0 



10:46:8.0 
R.A. {J200a) 



10:46:6.0 



Fig. 2. — Magnified image of tfie neutron star field 
{blue: a Chandra ~K.-v ay image after the PS'_F image 
reconstruction, red: a deep H band image taken 
witfi the IRSF observatory). The position of the 
X-ray source EHG7 is marked with a green circle 
with a 1" radius; it has no infrared counterpart. 
Another X-ray source about 10" southwest from 
EHG7 has a counterpart in the _ff-band. 



We then converted the source counts to energy flux 
using the PIMMS tool, assuming blackbody radia- 
tion with kT -128 eV and TVh -3.2x1021 cm'^, 
derived from a best-fit model to the XMM-Newton 
spectrum as described below. For the other data 
sets, we measured source flux from time averaged 
spectra, assuming the same spectral model. The 
results are shown in Table [1] EHG7 did not show 
any significant flux variation above —30% between 
the observations though a formal fit of these data 
rejected a constant flux model at above 90% con- 
fidence (xVd.o.f = 38.45/23). The X-ray emis- 
sion was stable for —30 years with i^x (0.3—2 keV) 
«1.06xl0~i'^ ergs cm^^ s^-^. 

Individual observations did not have sufficient 
photon statistics for detailed spectral analysis. 
Since the spectral shape did not change signifi- 
cantly between observationf@, we combined all the 



^Wc simultaneously fitted 0.3—3 keV spectra of all the avail- 
able EPIC cameras in each observation by an absorbed 
blackbody model and found that 95% confidence ranges 
of the Wh and kT parameters overlapped with each other. 



source spectra taken between 2000 and 2006 with 
the XMM-Newton MOS and pn detectors (Fig- 
ure[3]) . We used the HEAsoft tool mathpha to com- 
bine the source and background spectra, marf rmf 
to merge arf and rmf files and then addrmf to 
average produced response files over by weight- 
ing them by net source counts. The combined 
X-ray spectrum showed significant emission only 
below 2 keV and is almost featureless except for 
a slight dip at around 0.9 keV and a strong dip 
at around 0.6 keV. The latter feature at —0.6 keV 
was perhaps produced by edge absorption of oxy- 
gen in the interstellar medium and in the CCD 
detector response. The spectrum can be repro- 
duced using a 1-temperature optically-thin ther- 
mal plasma emission model suffering absorption 
by neutral gas along the line of sight (Table [2]), 
but the best-fit model has an unusually small el- 
emental abundance of <10~^ solar, without evi- 
dence of emission lines from oxygen, iron, or neon 
atoms near 0.5—1 keV. The spectrum also can be 
fit by a 1-temperature blackbody emission model 
with a temperature of 128±7 eV (90% confidence) 
suffering absorption by neutral gas with hydrogen 
column density of 3.2±0.4x 10^^ cm'^ (90% confi- 
dence). The spectrum could not be adequately fit 
using an ab s orbed power law model. We note that 
IPires et"aL ( 20081 ) derived a similar best-fit result 
in individual fits of XMM-Newton spectra. 

The hydrogen column density is consistent 
with interstellar absorption to t he Carina Nebula 
(— 3x 10^1 cm~2, see discussion in lLeutenegger et al. 



20031 ). The 90% confidence upper limit (3.6x lO^i cm'^) 
is much smaller than ab sorption through ou r 



Kalberla et al 



m ou r 
l2QQ5|), 



Galaxy (-1.2x10^^ cm 
which emission from an AGN should suffer. It is 
also significantly smaller than absorption to the 
background stars in the field, which have visual 
extin ctions twice as large as Carina cluster mem- 
bers (jPeGioia-Eastwood et al.ll200lh . This result 
suggests that EIIG7 is in front of the molecular 
cloud that lies behind the Trumpler 14 and 16 
clusters, and neighbors the Carina Nebula. 

The ratio of the X-ray flux to the bolomet- 
ric flux, log -Fx/^boi > —0.7, was significantly 
larger than that typically seen in stellar X-ray 
emission (< -3. iGagne et al.lll995l ). This large 
Fx/i^boi means that EIIG7 is unlikely to be a star 
in the Carina star-forming complex nor in the fore- 
and background but must be a compact object 
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Table 2 
Spectral Fits 



Model 


kT 




Abundance 






Ax' (d.o.f) 




(cV) 




(solar) 


(lO'^cm^^) 


(lO'^'crgs s"^) 




IT APEC 


171 (162- 


174) 


<8.0xl0— ' 


5.4 (5.0-5.8) 


32 


1.51 (125) 


IT Blackbody 


128 (121- 


133) 




3.2 (2.8-3.6) 


4.6 


1.49 (126) 



Note. — Absorption corrected Lx between 0.3 — 8 keV assuming d —2.3 kpc. Tlie parentheses in the 
kT and A^h columns denote the 90% confidence range. 



such as a black hole, a white dwarf or a neutron 
star. The X-ray characteristics of this source, a 
single temperature blackbody spectrum without 
time variation longer than hourly timescales, are 
typical of isolated neutron stars, where the emis- 
sion comes from the cooling neutron star surface. 
It is unlikely to be a black hole, which would 
show strong time variation on both short and long 
timescales, reflecting activity around the accre- 
tion disk near the event horizon. It is also un- 
likely to be a white dwarf, which has normally 
significa ntly lower blackbody temperature s than 
100 eV (jKahabka fc van den Heuvell Il997t ). Al- 
though some white dwarfs with masses close to 
the Chandrasekhar limit can have kT >100 eV, 
they should be more than five orders of m agni- 
tudes more luminous ( Sala fc Hernand 2005[ ) . All 
this evidence shows that the source is most likely 
a neutron star. 

We found no radio counterpart of EHG7 
in a continuum radio map of the Carina re- 
gion from the Southern Gal actic Plane Survey 
(|McClure-GrifEths et ani2005h . but the resolution 
of these data is only 1' and weak point sources 
will not be detected against the strong nebular 
radio sourc es Car I and Car II The ATNF Pulsar 
Catalogue ( Manchester et al. 20051 ) lists 11 pul- 
sars within 2° of EHG7, but no source coincident 
with EIIG7 itself. We searched for a pulse from 
the XMM-Newton pn data between 0.3-2 kcV 
in XMM060131 using the Fast Fourier Transform 
and searches. However, we detected no signif- 
icant pulses from EHG7 in the X-ray data in the 
frequency range between lO^-^-lS.e Hz. A 90% 
confidence upper limit on any signal amplitude in 
this frequency range was 27.7% (rms). 



s - 

1 X 



+ lwt|* 











1 2 
Energy (keV) 



Fig. 3.— Top panel: XMM-Newton MOS (black) 
and pn (red) spectra combined from all the avail- 
able data. The solid line shows the best-fit model 
of the spectra between 0.3—3 keV by an absorbed 
blackbody model. Bottom panel : the residuals of 
the fit. 

4. Discussion 

The intrinsic X-ray luminosity is 4.6 x 10'^^ ergs s"-'^ 
or equivalent to an emitting radius of ^^7.3 km 
assuming kT =128 eV and d =2.3 kpc. Assum- 
ing the stand ard cooling cu rve of a M =1.4 Mq 
neutron star I Tsuruta 20061 ). the age of EHG7 is 



estimated at ~0.5— 1x10 -year. The blackbody 
temperature, 0.13 keV, is also typical of neutron 
stars older than 1000 years (e.g. 1 E1207.4-5209 
PSR B0656-f 14 and PSR B 1055-52, ISanwal et al" 
2005f) . The X-ray spectrum 



2002; De Luca et al 



did not show a non-thermal power-law component 
in the hard band that would typically originate 
from charged particles in the magnetosphere. The 
upper limit of Lx [0.5-10 keV] < 5x10^° ergs s^^ 
estimated from a fit of the spectrum above 2 keV 
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assuming an additional power-law component with 
a fixed F =1.7 is similarly lower than is typical 
for middle-age d rotation-powered pulsars of a few 
xlO^ year old ( De Luca et al.|[2005l : lManzah et all 
20071 ) . These results suggest that the neutron star 
originated in a supernova explosion about 10^ 
years ago. 

The dip feature at ^--^0.9 keV can be produced 
by fundamental cyclotron absorption by elec- 
trons in the magnetic field with ^7.5x10^" Gauss, 
which is t ypical of middle-aged radio pulsars (e.g. 
Reisenegger..2001i) . or by protons in the magnetic 
fields with -^10^^ Gauss, reminiscent of the X- 
ray Dim Isol a ted Neutron Star (XDINS, see e.g. 
Haberll l2007l : [Kaplan ,200i). These stars have 
short pulsation periods of ~a few hundred milli- 
second and/or small pulse fraction of <20%. The 
existing X-ray data sets do not have enough time 
resolution nor photon statistics to detect such a 
pulse. 

Two other neutron stars hav e been discovered 
near other star fo rming regions ( Figer et al. 2005t 



of star formation. This is consistent with stel- 



Muno et al. 2006f l. Both of them have character- 
istics of magnetars — neutron stars with strong 
magnetic fields up to 10^^ Gauss — whose pro- 
genitor exploded within the last ~10* years and 
who might represent the high mass end of stars 
born in these star forming regions. However, the 
neutron star discovered in the Carina Nebula can- 
not be co-eval with the current generation of stars. 
Since the higher mass stars evolve faster, the pro- 
genitor of the neutron star would have to be more 
massive than r] Car. Eta Car's e stimated initial 
mass >150 Mq (jHillier et al.ll200lh is much higher 
than the conventional progenitor mass of neutron 
stars (less than ~25 Mq) and explosion of such 
a massive star would likely produc e a black hole, 



i nstea d of a neutron star although IWooslev et al 



(j2002l ) found that explosions of stars with masses 
up to ~60 Mq might produce a neutron star under 
certain circumstances. 

Based on standard stellar evolution models, the 
progenitor of the neutron star should have been 
born 6—30x10^ year ago with an initial mass 
between 8—25 Mq. This is significantly earlier 
than formation of the massive stars in the Trum- 



3x10^ years ago ( 


DeGioia-Eastwood et al. 2001; 


Massev et al. 2001 


). This suggests that the Ca- 



lar p opulation studies (jPeGioia-Eastwood et al 



20011 ) which also suggest that intermediate-mass 
stars have formed continuously over the last 10^ 
years. An alternative scenario is that EHG7 es- 
caped from a nearby star forming region, coming 
close to the Carina nebula by chance. One possible 
site is the Carina flare region located at ~2— 10° 
above the Carina nebula in Galactic coordinates. 
This region is suspected to have experi enced major 



star forming activity over 2x10^ yr (jFukui et al 
1999f ). Another possibility is that EHG7 ther- 



malizes through Bond i-Hoyle-Lyttleton accr etion 
from the nebular gas (jBlaes fc Madaul[T993l ). In 
that case, it can be much older than 10^ yr and 
its birth place is unpredictable. 

The last episode of star formation in the Carina 
Nebula appears to have been nearly contempora- 
neous with the supernova explosion of the progeni- 
tor of this neutron star. This may suggest that the 
expanding H II region, wind-blown bubble and/or 
supernova explosion of this star (and perhaps oth- 
ers yet undetected) played a role in triggering the 
last episode of star formation. Although the neu- 
tron star is currently well away from the nebula 
center, it could have moved the projected distance 
in ^10^ year with a transverse velocity of only 
~6 km s~^, which is much lower than typical kick 
velocities of ra dio pulsars (several hundred km s~^, 
Lvne fc Lorim cr 199^). 

Discovery of a neutron star supports the ar- 
gument that the diffuse high energy emission ob- 
served in the Carina Nebula originated in an an- 
cient supernova explosion. Some measurements 
require multiple supernova explosions in this field 
( Yonekura et al.ll2005t iHamaguchi et al.ll2007l ). so 



that the neutron star may be one of many neu- 
tron stars and black holes hidden in the Carina 
Nebula. If the soft extended X-ray plasma in real- 
ity heated up by supernova explosions ~10^ years 
ago, it would represent an important phase of the 
evolution of inter stellar medium when hot gas 
produced by supernova remnants merge together 
to form the hot ionized intercloud m edium (HIM) 
component (jMcKee fc Ostrikeiill977[ ). This result 
may have implications for the origin of diffuse X- 
ray emis sion observed now from many star f orming 
regions (|Townslev et~aLll2003HGiidel et al.ll2b08i) . 



rina Nebula has experienced at least 2 episodes 
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